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We  have  developed  an  approach  for  numerical  Monte  Carlo  solution  of  the 
Wigner  transport  equation  for  a  mesoscopic  system.  In  particular  the  Wigner j 
function  formalism  has  been  used  to  analyze  the  quantum  dynamics  of  charge 
carriers  in  presence  of  electron-phonon  interaction  in  open  mesoscopic 
systems,  making  use  of  the  perturbation  theory  with  the  external  fields 
incorporated  into  the  unperturbed  Hamiltonian.  Wigner  trajectories  have 
been  identified  both  for  the  ballistic  propagation  and  for  the  case  in  whict 
electron-phonon  scattering  is  included.  An  equation  for  the  Wigner  function] 
in  a  finite  region  inside  given  boundaries,  using  suitable  boundaries 
conditions,  has  been  derived.  As  a  case  of  study  a  step  potential  profile 
has  been  considered.  The  current  as  a  function  of  the  step  potential  as 
obtained  with  the  Wigner  function  corrected  by  the  effect  of  an  electron- 
phonon  scattering  process  is  presented.  On  the  coherent  transport  standpoidt, 
the  role  of  the  magnetic  field  on  the  condactance  of  wells  has  been  investijgated 
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Abstract 


We  have  developed  an  approach  for  numerical  Monte  Carlo  solution  of  the 
Wigner  transport  equation  for  a  mesoscopic  system.  In  particular,  the  Wigner 
function  (WF)  formalism  has  been  used  to  analyze  the  quantum  dynamics  of 
charge  carriers  in  presence  of  electron-phonon  interaction  in  open  mesoscopic 
systems  making  use  of  the  perturbation  theory  with  the  external  fields  in¬ 
corporated  into  the  unperturbed  Hamiltonian.  Wigner  trajectories  have  been 
identified  both  for  the  ballistic  propagation  and  for  the  case  in  which  electron- 
phonon  scattering  is  included.  The  concept  of  Wigner  trajectories  has  been 
used  as  a  guiding  line  to  face  the  problem  of  finding  an  equation  for  the  WF  in 
a  finite  region  inside  given  boundaries,  using  suitable  boundaries  conditions. 
As  a  case  of  study  a  step  potential  profile  has  been  considered.  The  current 
as  a  function  of  the  step  potential  as  obtained  with  the  WF  corrected  by  the 
effect  of  an  electron-phonon  scattering  process  is  presented.  On  the  coherent 
transport  standpoint,  the  role  of  the  magnetic  field  on  the  condactance  of 
wells  has  been  investigated. 
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I.  INTRODUCTION 


A  rigorous  quantum  transport  theory  of  electrons  in  semiconductor  structures  including 
both  coherent  propagation  and  scattering  mechanisms  is  still  lacking  even  though  its  need 
has  been  recognized  since  many  years.  In  particular  we  miss  a  theory  feasible  for  numerical 
calculations  to  be  compared  with  experimental  results  and/or  with  results  of  the  semiclassical 
theory  with  the  purpose  of  enlightening  the  different  predictions  of  the  two  models. 

The  Wigner  function  (WF)  approach  [1]  provides  a  rigorous  quantum-mechanical  tool 
and  constitues  a  direct  link  between  quantum  and  classical  descriptions  of  the  evolution  of 
the  system  in  phase  space.  In  this  approach  however  electron-phonon  (e-p)  interaction  has 
been  treated  using  either  relaxation  time  [2,3]  or  a  classical  Boltzmann  collision  operator 
[4].  In  this  work  several  advancements  in  the  efford  to  solve  the  above  problem  exactly  are 
presented. 

The  definition  of  the  WF  can  be  extended  to  include  electrons  interacting  with  phonons 
as  follows  [5]: 

fw  {r,p,ng,n'g,t)  =  j  (r  +  r72,  n,;  r  -  r72,  n'J  (1) 

where  r  and  p  are  the  electron  position  and  momentum,  p  is  the  density  matrix  for  the 
system  formed  by  an  electron  and  the  phonon  gas,  and  Uq  is  the  occupation  number  of  the 
phonon  mode  q.  Trace  over  the  phonon  variables  will  lead  to  the  traditional  electron  WF. 

In  what  follows  we  discuss  an  integral  equation  for  the  above  WF  and  apply  it  to  study 
quantum  electron  transport  in  open  systems  in  presence  of  scattering.  In  particular  suitable 
boundary  conditions  can  be  used  which  restrict  the  analysis  of  the  WF  inside  a  finite  domain. 
Wigner  trajectories  have  been  identified  also  for  the  case  of  e-p  interaction  and  constitute 
a  guiding  criterion  for  replacing  initial  conditions  with  boundary  conditions  for  the  integral 
equation. 

Our  theoretical  model  allows  us  to  evaluate  the  correction  to  the  WF  due  to  the  effect 
of  an  e-p  scattering  event. 
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II.  THEORETICAL  MODEL 


A  three-dimensional  system  of  indipendent  electrons  interacting  with  phonons  is  consid¬ 
ered  with  translational  invariance  along  two  directions  {x-y).  The  unperturbed  Hamiltonian 
Ho  of  the  system  contains  the  electron  Hamiltonian,  including  the  potential  profile  and  the 
free-phonon  term: 

Ho  =  He  -I-  Hp  =  ^  -h  V  (2)  -f-  ^  (2) 

where  p  is  the  electron  momentum,  m  is  the  electron  effective  mass,  V  (2)  is  the  electron 
potential  profile  (including  the  applied  voltage),  a,  and  af  are  the  annihilation  and  creation 
operators  of  the  phonon  q  with  frequency  ujq.  Hg.p  is  the  e-p  interaction: 

He_p  =  ^  ihF  (g)  -  a^ (3) 

Q 

where  F  (q)  is  a  function  that  depends  on  the  type  of  e-p  interaction.  In  the  present  work 
polar  optical  phonons  have  been  considered. 

The  staring  point  of  our  theoretical  approach  is  the  WF  generalized  to  include  phonons 
in  the  system,  in  a  single  electron  approach,  as  described  in  Eq.(l). 

The  time  evolution  of  the  density  matrix  is  described  by  the  Liouville-von  Neumann 

equation.  Starting  from  this  equation  we  derived  the  following  integral  equation  for  the 
WF: 


fw  {r,  p,  n„  n'^,  i)  =  ^  E  /-'  (r,  p)  f  dr'  [  dp'  {r',  p') 

n,n'  J 

xfw  {r',p',  rig,  n'^,  -|-  E  ^  (r,p)  dt' 

^  n,n' 

X  I  dr'  J  dp'  {n'  (nrig,  mm,)  {r',p')  (r',  p',  m„  n^,  t') 

TThTTlfj  ' 

~fnm  (^^  P')  fw  {r',  p' ,  n„  m„  t')  %'  {mmg,  n'nQ  |  (4) 


In  the  above  equation  'H'  —  and  H'  {nng,mmg)  are  its  matrix  elements  on  the 

basis  {\(l>i,ng)}  of  Ho.  The  coefficients  fim  (r,  p)  are  given  by 
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(5) 


fim  ir,p)  =  I dr'e-’f”-'  (r  +  r'/2\ 4>i)  |r  -  r' / 2) 

They  allow  to  move  from  the  WF  to  the  density-matrix  representation  and  viceversa.  Finally 
hu}  (n,  Uq)  indicates  the  total  unperturbed  energy  of  the  e-p  system  when  the  electron  is  in 
the  n-th  eigenstate  and  the  phonon  bath  is  in  a  state  with  n  phonons  in  mode  q. 

In  Eq.  (4)  the  first  term  in  the  r.h.s.  describes  the  ballistic  coherent  propagation  from 
the  initial  time  to  to  the  observation  time  t.  The  second  term  describes  the  contribution  of 
the  unknown  WF  to  the  “last”  interaction  vertex  at  any  time  t'  between  to  and  t  followed  by 
the  free  propagation  from  t'  to  t.  Even  though  the  equation  is  written  for  the  entire  system 
and  contains  also  the  phonon  variables,  it  has  been  shown  that  a  Monte  Carlo  solution 
allows  for  a  proper  average  over  the  phonon  coordinates,  as  long  as  hot-phonon  effects  are 
neglected  [6]. 


III.  WIGNER  TRAJECTORIES 

It  is  well  known  that  a  point-like  electronic  WF  in  phase  space  cannot  represent  a  real 
physical  system,  since  it  violates  the  uncertainty  principle.  We  may  however  consider  the 
evolution  in  time  of  a  (5— like  contribution  to  the  WF  5^  (r  —  Tq)  {p  —  p^)  fp  (uq,  nQ  inside 
Eq.  (4).  For  free  electrons  its  trajectory  in  phase  space  is  the  same  as  for  semiclassical 
particles.  In  fact,  for  a  plane-wave  basis  the  coefficients  fim  are  given  by 

(r,  P)  =  (6) 

and  the  ballistic  evolution  of  the  WF  as  given  by  the  first  term  on  the  rhs  of  Eq.  (4)  yields 

[7]: 

fw  (r,  p,  Uq,  n'q,  tj  =  fw  (r  -  ^(t  -  to)  ,p,  Uq,  n'g,  to'^  (7) 

where  m  is  the  electron  effective  mass  and  hu  [uq)  is  the  total  energy  of  the  phonon  bath 
in  the  state  {n,}.  The  first  factor  and  the  exponential  factor  describe  the  free  trajectory  of 
the  electron  and  the  time  evolution  of  the  free-phonon  bath,  respectively. 
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This  result  is  not  surprising  since  in  absence  of  scattering  and  up  to  quadratic  potentials 
(harmonic  oscillator)  a  differential  equation  can  be  written  for  the  WF  that  coincides  with  the 
Boltzmann  equation,  so  that  the  dynamical  evolution  of  the  WF  follows  the  same  trajectories 
of  the  representative  points  of  a  classical  gas.  This  implies  a  deformation  of  the  WF  while 
evolving  in  time  due  to  the  fact  that,  as  time  increases,  the  higher-momentum  components 
move  faster  than  the  lower-momentum  components. 

If  phonon  scattering  is  included,  similar  trajectories  exist,  with  very  interesting  proper¬ 
ties,  as  long  as  a  single  interaction  diagram  is  considered  without  time  integrations.  In  order 
to  prove  this  result  let  us  calculate  first  the  matrix  elements  of  the  interaction  Hamiltonian 
%'  in  Eq.  (4).  Assuming  a  constant  potential  and  using  again  plane  waves  as  basis  set,  we 
have: 


%'  {nng,mmg)=  {(l>n,ng\'H'  \(f>m,Tng) 

=  53  53  ^  (9)  n  V  ^  (1  -  {K  -  Oq  -  km)  ,  (8) 

(7=1, -1  q  q'jtg  V  ^ 

where  F  (g)  is  the  e-p  coupling  function.  Here  the  products  are  taken  over  all  values  of  q' 
except  q'  —  q. 

Substituting  Eqs.  (6)  and  Eq.(8)  in  Eq.(4)  and  using  the  (5-function  in  Eq.(8)  to  perform 
integration  over  the  electronic  wavevector,  we  obtain  for  the  second  term  in  Eq.(4): 

A/.  {r.p,n„n'„t)  =  ^  (p  -  f  df 

(27r)  „  V  l  j  Jto 

-i(a((n,)-u((n',))(t-t')  ^  /  dr'  f  dp'  ^  ^  E  (g) 

^  J  _ 1  1  /T. 


xe 


cr=l,— 1  q 


TI/oTTIq  — (T 


q'^q 

^^-iikn-<Tq-K,yr' ^3  (p'  _ 


''  llqlliq 

f  *  (^n 


fw  [r',p',mg,n'^,t') 


_e-(fen-fc„(-<7,)V^3  _  ^(fen  +  fcn'+^g)^  (r',p',ng,mg,t') 


f^q  d"  (1  T  (j)  Sm^in'^iSn'^mq+cr 


(9) 
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Performing  integrations  over  electron  coordinates  and  summing  over  the  electronic  states 


we  get: 


AU(r,p,n„n’^,t)  =  ^  ^ 

•/to  vn _ 1  1  ^ 


rriq  (T=ly-l  q 


1 


naTna—cr 


xfw(r-  ^{t  -  t')  ,p  -  aY,^q,n'g,t^ 

-fw  (r  -  ^{t-t'),p  +  ay, n„ m„ 


TTlq 

+  2  ^n'gmq+a  /  • 


(10) 


Substituting  Eq.(lO)  into  itself  we  obtain  the  second-order  term  corresponding  to  one  scat¬ 
tering  event.  The  second-order  correction,  e.g.  for  the  terms  corresponding  to  phonon 
emission,  results  to  be 


AU  {z,  p„  n„  n„  t)  =  2Sft{  ^  (q)  (n,)  dt'  ^  df 

x/»  t')  -  {t'  -  t")  -  ^  {f  -  to)  ,p,  +  hq,  n,  -  1,  n,  -  1,  Q  }  (11) 


The  trajectories  described  by  the  above  equation  are  represented  in  the  diagrams  con¬ 
tained  in  Fig.  1.  They  correspond  to  the  semiclassical  trajectories  where  q/2  is  transferred 
to  the  electron  at  each  vertex.  As  a  result  of  the  time  integrations,  however,  a  sum  over 
infinite  trajectories  gives  the  WF  correction  corresponding  to  a  single  scattering  event. 


IV.  BOUNDARY  CONDITIONS 

Guided  by  the  idea  of  Wigner  trajectory  introduced  in  the  previous  section,  we  face  the 
problem  of  finding  an  equation  equivalent  to  Eq.  (4)  for  the  WF  in  a  finite  region  inside 
given  boundaries.  In  fact,  to  solve  the  integral  equation  (4),  we  need  to  know  the  Wigner 
function  at  any  given  (2,^2)  at  time  to.  Such  initial  condition  is  not  suitable  for  describing 
an  open  system.  In  fact,  by  definition,  a  system  is  open  when  it  can  exchanges  particles 
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with  the  environment:  this  means  that  there  is  a  boundary  between  the  system  of  interest 
and  the  rest  of  the  universe.  Therefore  to  model  such  a  system  the  Wigner  function  has  to 
be  known  at  the  initial  time  to  inside  the  system  and  at  all  times  on  the  boundary.  Prom  the 
analytical  point  of  view  this  request  corresponds  to  substituting  in  Eq.  (4)  the  integral  over 
the  whole  space  by  three  terms.  One  term  is  identical  to  the  original  expression  in  Eq.  (4), 
with  the  integral  over  the  initial  condition  limited  to  the  points  internal  to  the  device.  The 
other  two  terms  are  evaluated  at  the  two  boundaries  of  the  device.  In  the  following  we  will 
derive  the  expression  of  this  terms  for  the  case  in  which  the  eigenstates  of  the  unperturbed 
Hamiltonian  are  plain  waves,  considering,  for  the  sake  of  simplicity,  the  one  dimensional 
case.  Under  this  assumption  the  coefficients  of  Eq.  (5)  can  be  rewritten  as 

fnn'{z,Pz)  =  J 

=  (p,  -  .  (12) 

Let  us  consider  a  device  defined  between  z  =  -A  and  2  =  +A.  We  will  analyse  now  the 
ballistic  term 


{z,P.,n^,n’^,t)  f  dz'  f  dp'J*^,  {z',p',) 

nn' 

•  (13) 

Substituting  Eq.  (12)  in  Eq.  (13),  considering  at  first  -oo  <  2'  <  -A,  leads  to 

fw^  Pz^ 

=  2  [  f  dkn'  (pz  —  ^  ^  {t~to  )  Xlq  (”9-^q  )^q  (<“*0  ) 

h  (27r)  J  J  y  2  y 

I  (p'^  -  u  - 

y- 1  dz"f  /.  . 

From  the  properties  of  the  ci-function  we  get 


(14) 
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(15) 


z'  <-A=^z-^{t-to)<-A 
m 


t>  to  +  P'z  >  0 

m 


Since  in  this  specific  case,  particle  entering  the  device  from  the  left  boundary,  the  condition 
<  0  can  not  be  verified,  we  are  left  with 


^ 

to<t - i  <  t,  if  p;  >  0. 


El 

m 


(16) 


Setting 


~  ^  d'  —  ^  - 


m 


m 


m 


(17) 


yields 


z' 


m 


(18) 


Then  we  can  perform  a  change  of  coordinates  from  {z',p'^)  to  (r,  ^)  as  follows; 


< 

\ 

-x4  -  ^  (r  -  to),  to<T  <t 

< 

5 

det  J  =  det 

( 

1 

-  to) 

P'z=^  ,  ^>0 

K 

0  1 

1 

/ 

m 


(19) 


obtaining  in  conclusion: 

ft  .  /•+0O  _ 


Jto  JO  m 

^  w 

X  I dz"j ip"j;Az",pl)f,[z’,p'',n„n',,u) 

=  rEAn-l^.P^) 

^  nn'  ■’to 

X  ^/nn'  (-^1 P^)  /«.  p'z,  n^,  n'g,  If)  .  (20) 
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The  same  procedure  can  be  applied  in  the  case  +A  <  z'  <  +oo. 

The  details  of  the  derivation  for  the  terms  including  scattering  events  require  a  more 
sophisticated  analysis  and  are  discussed  in  Ref.  [7].  Here  we  simply  report  the  final  result: 


fw  (z,p^,ng,n'g,t)  = 

I E  (z,P.)  £  dz'  J  dp'X..  /.  {z'.  P'., 

E  (z.  p,)  £  dt'  { [  £”  dp',  i/;,,  (-.4,^;) 

xfw  {-Ap'z,  ng,  n'g,  t')  -  f  dp'^  {A,p'^)  U  (A,p'^,  n„  n'  t') 

+  dz'  J  dp[  fmn'  (2',  P'z)  fw  {z',p'z,  Tl'g,  t') 

f nm  ^Pz)  fw  ^Pzi  ^95  ^  "^q^  j  ^  * 


(21) 


where  again  for  simplicity  a  one-dimensional  case  has  been  considered,  and  —A  and  A  are 
the  device  boundaries. 

In  order  to  verify  the  correctness  of  the  use  of  the  boundary  conditions  for  an  open 
mesoscopic  system  we  have  first  evaluated  the  ballistic  evolution  of  the  first  term  in  Eq. 
(21)  assuming  a  step  potential  profile,  =  0  at  t  =  to  inside  the  device,  and  a  Maxwellian 
local-equilibrium  incoming  distribution  at  the  left  and  right  boundaries.  Results  are  shown 
in  Fig.  2  at  different  times  after  the  initial  condition. 

At  the  longest  times  we  recover  for  /„  the  same  result  which  is  obtained  without  using 
boundary  conditions  for  a  stationary  infinite  system,  where  /«,  is  evaluated  by  means  of  an 
equilibrium  density  matrix  diagonal  over  the  scattering  states. 


V,  APPLICATIONS 

In  order  to  evaluate  the  effect  of  an  e-p  scattering  event  on  the  WF  we  have  performed  an 
iterative  expansion  of  Eq.  (21)  up  to  the  second  order  in  the  e-p  perturbation  Hamiltonian. 
As  a  case  of  study  we  have  considered  a  step  potential  profile.  Integrations  in  time  have 
been  performed  analytically;  the  integrals  over  z',  p'^  and  over  the  electron  states  have  been 
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evaluated  with  full  numerical  integration  while  the  integral  over  the  phonon  modes  has  been 
performed  by  means  of  a  Monte  Carlo  selection.  Space  correlations  for  the  evaluation  of  the 
coefficients  in  Eq.  (5)  have  been  considered  up  to  40  nm.  The  steady-state  WF  for  a  0.1  eV 
step  including  the  effect  of  an  e-p  scattering  event  is  presented  in  Fig.  3.  Fig.  4  shows  the 
current  as  a  function  of  the  step  potential  as  obtained  with  the  ballistic  WF  (solid  curve)  and 
with  the  WF  corrected  by  the  effect  of  an  e-p  scattering  process  (circles)  switched  on  50  fs 
before  the  “observation  time” .  Comparison  is  presented  with  the  outcomes  of  a  semiclassical 
calculation  based  on  the  Boltzmann  equation.  The  quantum  ballistic  curve  is  always  lower 
then  the  corresponding  semiclassical  one,  and  instead  of  reaching  a  saturation  value,  rise 
up  to  a  maximum  and  then  decreases.  This  is  due  to  the  fact  that  in  the  quantum  picture 
increasing  the  potential  step  leads,  as  a  consequence  of  quantum  reflection,  to  a  decrease  of 
the  transmission  coefficient  that  is  to  a  decrease  of  the  current  intensity.  Finally,  in  both 
quantum  and  semiclassical  calculations  the  effect  of  the  e-p  interaction  is,  as  expected,  a 
reduction  of  the  current  with  respect  to  the  ballistic  case. 

VI.  COHERENT  TRANSPORT  IN  MESOSCOPIC  DEVICES 

Concerning  the  study  of  coherent  phenomena  in  open-boundary  2D  systems,  during  this 
last  year,  the  research  activity  has  proceeded  both  from  the  theoretical  and  the  application 
standpoints. 

As  for  the  first  one,  the  inclusion  of  the  magnetic  field  inside  the  numerical  solver  of  the 
2D  Schrodinger  equation  has  been  achieved.  This  has  allowed  to  investigate,  for  example, 
the  role  of  the  magnetic  field  on  the  conductance  of  rough  quantum  wells  [8].  While  the 
conductance  of  an  ideal  quantum  well  is  greately  deteriorated  by  the  roughness  present  at  its 
surface,  the  presence  of  a  magnetic  field  partly  recovers  it.  The  magnetic  field  modifies  the 
coherent  behaviour  of  electrons  leading  to  conduction  properties  that  depend  on  its  strength 
[9]. 

We  have  analized  a  rough  quantum  wire,  where  the  roughness  is  modelled  with  a  gaussian 
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autocorrelation  and  the  potential  is  zero  inside  the  channel  and  infinity  at  the  boundary. 

In  a  rough  quantum  wire  the  magnetic  field  reduce  the  effect  of  roughness  increasing  the 
conductance. 

Fig.  5a  shows  the  squared  modulus  of  the  electron  wave-function  for  the  case  of  an 
electron  in  a  rough  channel.  The  carrier  is  in  the  fundamental  subband,  since  it  is  supposed 
to  have  total  energy  lower  than  the  threshold  of  the  first-excited  subband,  and  is  coming  from 
the  left  side  of  the  channel.  Getting  close  to  the  right  boundary,  since  the  channel  becomes 
wider,  the  electron  moves  into  the  upper  subband.  This  produces  a  reflection  of  the  electron, 
and  decreases  the  conductance.  In  presence  of  a  magnetic  field  the  electron  approaches  one 
of  the  edge  of  the  channel  (Fig.  5b:  positive  magnetic  field;  Fig.  5c:  same  value  of  the 
field  as  in  Fig.  5b,  but  with  negative  sign)  and  the  backscattering  is  strongly  reduced.  The 
two  situations  shown  in  Fig.s  5b-5c  are  equivalent  from  the  transport  standpoint.  In  fact, 
because  of  the  Onsager  relation,  that  is  a  consequence  of  the  time  inversion  symmetry,  the 
conductance  depends  only  on  the  absolute  value  of  the  magnetic  field  not  on  its  sign. 

The  effect  of  the  magnetic  field  on  the  conduction  properties  of  rough  channel  is  indepen- 
det  from  the  roughness.  This  can  be  seen  studing  a  large  ensamble  of  rough  channels  with 
different  roughness.  In  Fig.  6  it  is  shown  the  distribution  of  the  conductance  values  over  the 
samples,  at  increasing  magnetic  field.  The  trasversal  region  in  wich  the  electron  is  costrained 
by  the  magnetic  field  is  well  indicated  by  the  magnetic  length  Ib  =  /fi/|eS|,  where  B  is  the 
value  of  the  magnetic  field.  The  label  L-inf  in  the  first  graph  of  Fig.  6  indicates  the  zero 
magnetic  field  case,  in  the  other  two  graphs  the  magnetic  field  increases,  and  the  magnetic 
length  Ib  takes  the  values  L/2  and  L/4,  where  L  is  the  mean  width  of  the  channel.  The 
conductance  distribution  at  zero  magnetic  field  shows  a  tail  at  low  conductance  values  that 
is  progressively  reduced  by  increasing  the  magnetic  field. 

On  the  application  side,  the  effect  of  channel  quantization  in  silicon  MOSFET  has  been 
investigated.  In  particular,  different  devices,  which  are  considered  possible  alternatives 
to  standard  MOSFETs  for  the  next  ULSI  generations,  have  been  explored.  In  [10],  [11] 
threshold  voltage  shift,  capacitance  reduction,  and  mobility  enhancement  have  been  studied 
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for  highly  non-homogeneous  devices.  In  these  devices,  a  low-doping  layer  is  epitaxially  grown 
on  a  standard  high-doping  substrate.  As  can  be  seen  in  [10],  [11],  the  presence  of  the  epi- 
layer  doesn’t  degrade  the  capacitance  of  the  MOS  while  mobility  is  improved  by  the  use  of 
low  dopings.  In  addition,  short  channel  effects  are  still  under  control  through  the  buried 
high  doping  layer. 

The  same  kind  of  investigation  has  also  been  performed  on  single-  and  double-gate  MOS 
devices  [12],  [13].  In  particular,  the  relevance  of  the  so  called  volume  inversion  due  to  size 
quantization,  as  opposed  to  the  reduction  of  transverse  effective  field,  in  improving  the 
electron  low-field  mobility  in  such  devices  has  been  explored. 


VII.  CONCLUSIONS 

In  this  final  report  we  have  presented  the  main  results  achieved  under  this  contract. 

A  theoretical  and  computational  analysis  of  the  quantum  dynamics  of  charge  carriers  in 
presence  of  electron-phonon  interaction  based  on  the  Wigner  function  has  been  applied  to 
the  study  of  transport  in  mesoscopic  systems. 

We  have  shown  the  existence  of  Wigner  trajectories  both  for  the  ballistic  evolution  of 
the  WF  and  for  the  case  in  which  e-p  scattering  is  included. 

We  have  faced  the  problem  of  finding  an  equation  for  the  WF  in  a  finite  region  inside 
given  boundaries.  We  have  proved  that  for  the  solution  of  the  integral  equation  describing 
the  WF  it  is  possible  to  substitute  the  knowledge  of  the  WF  over  all  space  at  a  given  (initial) 
time  to  with  the  knowledge  of  the  same  function  inside  the  region  of  interest  at  to  and  at  all 
times  at  the  boundary  for  “entering”  momenta. 

From  the  application  standpoint  we  have  analyzed  the  effect  of  an  e-p  scattering  event 
on  the  WF  for  the  case  of  a  step  potential  profile.  Results  for  the  current  in  comparison 
with  semiclassical  calculations  are  presented. 

Concerning  the  study  of  coherent  phenomena  in  open-boundary  2D  systems,  the  inclusion 
of  the  magnetic  field  inside  the  numerical  solver  of  the  2D  Schrodinger  equation  has  been 
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achieved.  This  has  allowed  to  investigate  the  role  of  the  magnetic  field  on  the  condactance 
of  rough  quantum  wells.  It  is  found  that  in  both  structures  the  magnetic  field  reduce  the 
effect  of  roughness  increasing  the  conductance. 

Finally  the  effect  of  channel  quantization  in  silicon  MOSFETs  has  been  investigated.  In 
particular,  different  devices  which  are  considered  possible  alternatives  to  standard  MOS¬ 
FETs  for  the  next  ULSI  generations,  have  been  explored. 

VIII.  PUBLICATIONS 

•  C.  Fiegna,  A.  Abramo,  “Optimization  of  channel  profiles  for  ultra-short  MOSFETs  by 
quantum  simulation”,  in  lEDM  Tech.  Digest,  815,  S.  Francisco,  CA  (USA),  Dec.  1996. 

•  A.  Abramo,  “A  general  purpose  2D  Schrodinger  solver  with  open/closed  boundary 
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Fig.l  Example  of  the  Wigner  trajectories  that  contribute  to  the  result  in  Eq.  (4)  for 
real  and  virtual  emissions  of  a  mode  q. 
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Fig.4  Current  as  a  function  of  the  step  potential.  Comparison  between  quantum  and 
semiclassical  calculations  with  and  without  e-p  interaction.  Wigner  ballistic  curve  (solid 
line),  Wigner  curve  with  e-p  scattering  (circles),  semiclassical  ballistic  curve  (dashed  line), 
semiclassical  curve  with  e-p  scattering  (dotted  line).  Carrier  density  at  the  contacts  n  =  10^® 
cm“^. 
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Fig. 5  Squared  modulus  of  tlic'  electTon  wave-fuiietion  in  a  rough  channel  for  the  case  of 
a  carrier  in  the  fuudaiiient.al  sul)band  having  total  energy  lower  then  the  threshold  of  the 
first  excited  subband  and  coming  from  the  left,  a)  No  magnetic  field;  b)  positive  magnetic 
field;  c)  nagative  magnetic  fic'ld  with  the  same  value  as  in  b). 
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Sample  distribution 
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Fig.6  Distribution  of  the  conductance  values  over  an  ensamble  of  samples,  at  increasing 
magnetic  fields,  normalized  to  the  peak  value  (see  the  text  for  the  meaning  of  the  symbols). 
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